REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  NO.  0704-0188 

Public  Reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comment  regarding  this  burden  estimates  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188,)  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  (  Leave  Blank) 

2.  REPORT  DATE  02/22/01 

3.  REPORT  TYPE  AND  DATES  COVERED 

Final  Project  Report:  9/96  -  9/00 

4.  TITLE  AND  SUBTITLE 

Modeling  the  Hydrolysis  of  Phosphonate  and  Phosphate  Esters 

5.  FUNDING  NUMBERS 

DAAH04-96-1  -0424 

6.  AUTHOR(S) 

Christopher  J.  Cramer 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Department  of  Chemistry  and  Supercomputer  Institute 

University  of  Minnesota 

207  Pleasant  St.  SE 

Minneapolis,  MN  55455 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

U.  S.  Army  Research  Office 

P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be  construed  as  an  official 
Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 

12  a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

12  b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

Final  summary  of  work  on  development  and  application  of  models  for  predicting  pathways  for  the  destruction  of  chemical 
weapons  agents  in  aqueous  solution. 

14.  SUBJECT  TERMS 

VX,  molecular  modeling,  detoxification,  solvation,  environmental  remediation 

15.  NUMBER  OF  PAGES 

6 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OR  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION 

ON  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

UL 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.2-89) 


Prescribed  by  ANSI  Std.  239-18 


298-102 


MODELING  THE  HYDROLYSIS  OF  PHOSPHONATE  AND  PHOSPHATE  ESTERS 


FINAL  PROGRESS  REPORT 

CHRISTOPHER  J.  CRAMER 

FEBRUARY  22,  2001 

U.S.  ARMY  RESEARCH  OFFICE 

DAAH04-96- 1-0424 

UNIVERSITY  OF  MINNESOTA 

APPROVED  FOR  PUBLIC  RELEASE; 

DISTRIBUTION  UNLIMITED 

THE  VIEWS,  OPINIONS,  AND/OR  FINDINGS  CONTAINED  IN  THIS  REPORT  ARE 
THOSE  OF  THE  AUTHOR(S)  AND  SHOULD  NOT  BE  CONSTRUED  AS  AN  OFFICIAL 
DEPARTMENT  OF  THE  ARMY  POSITION,  POLICY,  OR  DECISION,  UNLESS  SO 


DESIGNATED  BY  OTHER  DOCUMENTATION 


2 


1.  FOREWORD:  N/A 

2.  TABLE  OF  CONTENTS:  N/A 

3.  LIST  OF  APPENDICES,  ILLUSTRATIONS,  AND  TABLES:  N/A 

4.  SUMMARY  OF  RESEARCH 

A.  THE  PROBLEM  STUDIED:  We  proposed  to  develop  and  apply  new  computational 
chemistry  techniques  for  predicting  the  structure  and  reactivity  of  phosphorus-containing 
compounds  in  aqueous  solution.  Early  efforts  were  to  (i)  refine  a  basis  set  that  optimally 
balances  efficiency  and  flexibility  for  ab  initio  calculations  on  molecules  containing 
phosphorus,  (ii)  extend  our  quantum  mechanical  Generalized  Bom/Surface  Tension 
continuum  solvation  models  (the  SMr  series  of  models)  to  ab  initio  self-consistent  field 
implementations,  and  (iii)  begin  exploring  gas-phase  potential  energy  surfaces  for  hydrolysis 
and  nucleophilic  substitution  reactions  of  phosphonate  and  phosphate  esters  and  ester 
derivatives.  The  intent  of  the  latter  effort  was  to  provide  potentially  useful  starting  points  for 
calculations  including  the  effect  of  aqueous  solvation,  and  moreover  to  permit  an  evaluation 
of  exactly  how  important  those  effects  are  in  controlling  structure  and  reactivity. 

Midterm  focus  was  on  applying  the  completed  solvation  models  to  the  nucleophilic 
substitution  reactions  of  phosphonate  and  phosphate  esters  and  ester  derivatives.  Particular 
attention  was  paid  to  mechanistic  details  of  the  overall  reactions.  I.e.,  do  they  proceed  in  a 
concerted  fashion,  or  are  stable  intermediates  involved?  When  intermediates  are  involved, 
how  does  stereochemistry  associated  with  trigonal  bipyramidal  phosphoranes  affect  the  rate 
and  product  distribution  of  the  reaction?  How  do  phosphonates  compare  to  phosphates  in 
terms  of  reactivity  and  selectivity?  How  does  substitution  of  sulfur  for  oxygen  modify  the 
structure  and  reactivity  of  phosphono-  and  phosphorothioates? 

Final  efforts  built  upon  the  earlier  results  and  focused  specifically  on  aqueous 
hydrolysis  of  phosphonate  and  phosphate  esters  and  ester  derivatives.  Mechanistic  analysis 
required  a  careful  investigation  of  substrate  acidity  constants  and  activation  energies.  Specific 
mechanistic  details  associated  with  known  chemical  weapons  agents  were  addressed.  In 
particular,  attention  was  paid  to  the  individual  steps  in  perhydrolysis  of  VX.  The  relative 
tendencies  for  cleavage  of  different  P-X  bonds  in  known  nerve  agents  under  aqueous 
hydrolysis  conditions  were  characterized. 
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B.  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

During  the  period  of  this  contract,  we: 

•  examined  the  utility  of  pairwise  screening  algorithms  to  replace  complicated 
volume  integrals  with  sums  of  analytical  pairwise  terms.  This  approach  improves 
computational  efficiency  and  should  have  broad  application  to  molecular  simulations,  in 
addition  to  serving  with  our  Generalized  Bom  continuum  model. 

•  characterized  the  key  steps  in  the  peroxyhydrolysis  of  VX  nerve  agent  in 
aqueous  solution.  This  included  consideration  of  all  possible  intermediates  in  the 
addition/pseudorotation/elimination  pathway  at  suitable  levels  of  electronic  structure  theory 
including  solvation  effects  by  continuum  methods. 

•  developed  and  calibrated  a  theoretical  model  for  the  successful  prediction  of 
reaction  rates  for  hydrogen  atom  abstraction  in  aqueous  solution.  The  motivating  factor  was 
the  possible  utility  of  hydroxyl-radical-based  schemes  for  detoxification  of  chemical 
weapons  (e.g.,  from  ionization  of  liquid  water  using  electron-beam  technology). 

•  developed  improved  models  for  representing  solute  charges  in  aqueous 
solution.  Such  models  were  particularly  helpful  in  extending  our  calculations  into  the  ab 
initio  regime  so  as  to  permit  a  more  accurate  handling  of  chemical  weapons  agents 
containing  phosphorus. 

•  extended  our  continuum  solvation  models  to  the  Hartree-Fock  and  Density 
Functional  levels  of  theory.  At  these  ab  initio  levels,  we  demonstrated  much  improved 
accuracy  over  previous  semiempirical  levels  for  molecules  containing  second-row  atoms 
like  phosphorus. 

•  developed  and  coded  the  theory  for  analytical  gradients  including  continuum 
solvation  effects.  This  advance  permitted,  for  the  first  time,  optimization  of  structures  in 
solution,  greatly  facilitating  the  location  of  key  reaction  transition  states.  The  critical  issue  is 
that  such  transition  states  may  not  occur  in  the  gas  phase,  and  thus  they  cannot  be 
approximated  by  simply  solvating  the  gas-phase  structure — instead,  they  must  be  located  in 
solution. 


•  developed  and  continue  to  distribute  our  licensed  codes,  amsol  and  gamesol,  to 
calculate  solvation  free  energies  using  semiempirical  and  ab  initio  quantum  mechanical 
continuum  solvation  models.  The  codes  are  distributed  by  ourselves  and  various 
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commercial  vendors  under  license  from  the  Regents  of  the  University  of  Minnesota  (see 
technology  transfer  below). 

•  developed  the  first  quantum  mechanical  solvation  model  to  calculate  partitioning 
of  organic  molecules,  to  include  pesticides  analogous  to  chemical  weapons  agents,  between 
water  and  soil.  Such  models  are  critical  to  understanding  the  environmental  fates  of 
weapons  agents  insofar  as  environmental  detoxification  depends  on  the  environmental 
phase  in  which  agent  sequestration  takes  place. 

•  documented  the  utility  of  our  solvation  models  in  predicting  one-electron  redox 
potentials.  Again,  the  motivation  is  to  enable  modeling  of  alternative  (electrochemical) 
methodologies  for  the  destruction  of  dangerous  environmental  contaminants  like  chemical 
weapons  agents. 
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E.  TECHNOLOGY  TRANSFER 

We  developed  several  computer  programs  for  performing  semiempirical  and  ab 
initio  calculations  that  include  the  effects  of  aqueous  and  non-polar  solvation.  Details  are 
available  at  www.comp.umn.edu  on  amsol,  dgauss,  gamesol,  mn-gsm,  hondosol, 
zindosol,  and  omnisol.  Most  of  these  programs  are  all  available  for  free,  under  license,  to 
non-profit  researchers  and  are  distributed  by  the  University  of  Minnesota.  Various  vendors 
distribute  our  codes  to  for-profit  institutions  under  licensing  arrangements  with  the 
University. 
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